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Abstract 

We study Ginzburg-Landau equations for a complex vector order parameter ^ = 
(ip +,%/>-) € C 2 . We consider symmetric vortex solutions in the plane M 2 , ip(x) = 
f±(r)e m±d , with given degrees n± G Z, and prove existence, uniqueness, and asymp- 
totic behavior of solutions as r — > oo. We also consider the monotonicity properties 
of solutions, and exhibit parameter ranges in which both vortex profiles /+,/- are 
monotone, as well as parameter regimes where one component is non-monotone. The 
qualitative results are obtained by means of a sub- and supersolution construction and 
a comparison theorem for elliptic systems. 
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1 Introduction 

In this paper we study entire solutions in M 2 of coupled systems of two Ginzburg-Landau 
type equations, for complex vector-valued functions \I/ = (ip + (x), ip^(x)) : M. 2 — > C 2 , 

-A^p + + [A + (\^\ 2 - t\) + B{\^\ 2 - ti M + = 0,1 
-A^_ + [A_(\iIj_\ 2 - t 2 _) + B(\ip + \ 2 - tl)}tp_ = O.J 

Coupled systems of two Ginzburg-Landau equations of the type studied here arise both 
in the modeling of p-wave superconductors [K nRo] and in two-component Bose-Einstein 
condensates (BEC) |EKNTT] . We give a brief discussion of the connection between the 
energy ^ and the two component BEC model at the end of this section. 

Throughout the paper we make the following assumptions concerning the constants ap- 
pearing in Q: 

A + ,A_ > 0, B 2 < A + A„, t+,t_>0. (H) 
The system ([I]) gives the Euler-Lagrange equations formally associated to the energy 

= f ||V\I/| 2 + ^[^ + (|^ + | 2 -^) 2 + ^_(|^_| 2 -t 2 _) 2 + 2 J B(|^ + | 2 -^)(|^_| 2 -t 2 )] (2) 

for functions \I/ = ipJ) G H 1 {VL; C 2 ). As is the case for any Ginzburg-Landau model, this 
energy is divergent for Q = M. 2 for nontrivial entire solutions. Nevertheless, the problem is 
variational in nature, and the functional E(ty; Q) will be a useful tool for studying solutions 
locally. 

By hypothesis (H), the potential term in the energy 

F(¥) = A + (\i; + \ 2 - tl) 2 + A.{\^\ 2 - tlf + 2B(\iP + \ 2 - t 2 + )(|^| 2 - t 2 _) 

is strictly positive definite, and attains its minimum (of zero) when \^±\ = t±. Following 
previous work on the classical Ginzburg-Landau equations (see |BMR] ) . we seek solutions 
for which the potential energy is integrable, 

F(^) dx < oo. 

By the positive definiteness of F, this suggests that the solutions we seek should attain the 
asymptotic values 

|^ ± (a;)| —> t± as \x\ — > oo. 

That is, for \x\ large, the solutions may be written in polar form, ip±( x ) = p±(x)e l ^ >± ^ with 
real-valued p±(x), <f)±(x), and p±(x) ~ t±. The phases <j)±(x) may have nontrivial winding 
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number around any large circle Cr enclosing the origin, and thus these solutions carry two 
integer degrees n± = deg(-n^r , C#) G Z. As there are no smooth S-valued functions with 
nontrivial degrees in a simply connected domain, solutions with nontrivial winding 

must vanish in one or more of its components ip± to avoid singularities. These zeros are the 
vortices of the solution. 

In this paper we consider special solutions of the equations ([I]) with vortices, which are 
equivariant with assigned degree pair G Z 2 , 

in polar coordinate (r, #). By taking complex conjugates if necessary, we may assume that 
n± > 0. By the equivariant ansatz, the system ([!]) reduces to a system of ODEs: 
1 n 2 

2 

n_ 



(3) 



-/+ - "/+ + 7?/+ + [Mf+ - *+) + B{f- - £)] U = 0, for r G (0, oo), 
1 n 2 

-f- ~ -f- + -^f- + [A-U- ~ t-) + BUI ~ tl)] f- = 0, for r G (0, oo), 

/±(r) > for all r G [0, oo), 
/±(r) — > t± as r — > oo, 
/±(0)=0 if n ± ^0; /4(0) = if n± = 0. 
The first result is the existence of unique equivariant solutions, for any given degree pair 
[n+,n_] G Z 2 : 

Proposition 1.1. Lei n± G Z 6e ^iuen and A + A_ — B 2 > 0. T/ien £/iere exists a unique 
solution [/+(r), /-(?")] to ([3]) /or r G [0, oo). Moreover, 

f ± G ((0, oo)) , (4) 
/±(r) > /or a// r > 0, (5) 
/ ± (r) ~r |n±l /orr ~ 0, 

/■oo 

/ (f'±) 2 rdr < oo, 
./o 

/ [A + (/ 2 - t 2 ,) 2 + 2B(f 2 + - t 2 + )(f 2 - t 2 _) + A4f 2 - tl) 2 ] r dr 
Jo 

In particular, ty(x) = \f + (r)e m+e , f_(r)e m ~ e ] is an entire solution of Q in K 2 . 

We note that conclusion ^ gives a quantization result (see jBMRj ) for the solutions in 
terms of the degree pair at infinity, 



n 2 t 2 +n 2 _t 2 _. 



(6) 
(7) 

(8) 



I 2 - t 2 + y + a_(|^_i 2 - 1 2 ) 2 + 2s(i^ + i 2 - d)(iv-i 2 - t 2 ji 



2tt (n 2 t 2 +n 2 t 2 ) . (9) 
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To determine the shape of the vortex profiles, we first consider the asymptotic form of 
the solutions for r — > oo. We prove that 

/±(r)=t ± + ^! + ^ + 0(r- 6 ), /4(r) = -^ + 0(r- 5 ), as r ^ oo, (10) 

with 



a± 



1 Bn 2 T - A T n± 

2 {A + A_ - B 2 )t±' 



and (rather complicated) constants b± given in (27). A formal asymptotic expansion of this 
form (in fact, an expansion to arbitrary order in 1/r) may be obtained by simply substituting 
an ansatz into the system of equations and matching terms. Our results, presented in 



Theorem provide rigorous confirmation of this expansion by means of a sub- and super- 
solution construction. Moreover, we prove that the expansion is uniform in the coefficients 
A±, B, t± lying in a compact set. The proof is completed using a new and original comparison 



principle (Lemma 3.1) for elliptic systems, which generalizes the one in |AB]. While there is 
no comparison or maximum principle which applies to a generic elliptic system, in this case 
the equations have a special structure (similar to competitive or cooperative systems) which 
fortunately admits a comparison result. 



From the asymptotics (10), we see that the shape of the solutions depends strongly on 
the coefficients, in particular the sign of the interaction coefficient B. For B < 0, both 
components approach their limiting value t± from below, as is familiar from the classical GL 
vortices. However, for B > 0, this may no longer be the case, and for certain choices of n± 
and B one of the components will approach its limiting value from above. Such behavior was 
already noted in |ABMlj in the case n_ = 0, for a "balanced" system, A + = A_, t + = t_. 
In case n± ^ 0, our result implies that there are parameter regimes in which vortex profiles 
will be non-monotone. 

Finally, we consider the question of monotonicity of the solution profiles. For the standard 
GL vortices, the vortex profile is known to be strictly monotone increasing in r. As suggested 
by the asymptotic expansion above, the validity of this property is strongly dependent on 
the value of B. We prove the following: 

Theorem 1.2. Let A + ,A^ > be fixed, and B such that B 2 < A + A-. Assume ty(x;B) = 
[f+{r; B)e m+d , /_(r; B)e m - e ] is the equivariant solution for those parameters A±, B. 

(i) If B < 0, then f'±(r; B) > for all r > and for any degree [n+, n_] . 

(ii) If B > 0, n + > 1 and n_ = 0, then f' + (r; B) > and f'_(r; B) < for all r > 0. 
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(iii) For any pair [n + ,n_] with n + ^ ^ n_ ; i/iere exists B > swc/i i/wrf f ; ± (r;B) > 
/or all r > and all B with < B < B . 

As it is expected that only degrees n± = 0, ±1 give rise to locally minimizing solutions 
(see |Shj . |ABM2j ). let us summarize what we know about these specific symmetric vortex 
solutions. The situation for the [n+,n_] = [1,0] vortex is very clear: when B < both 
components are monotone increasing, but when B > 0, /_ decreases monotonically, while 
/+ increases. This is the same result as was found in |ABMlj in the special case A + = A_, 
t+ = t-. 

For the [n+,n_] = [1,1] vortex the picture is more complicated. For the interaction 
coefficient B < 0, both components f±(r) are monotone increasing to their limiting values 
t± as r — > oo, and so they always resemble classical GL vortices. When B = 0, the equations 
decouple, and each component is a rescaled GL vortex profile, and hence monotone increasing 
as well. For B > 0, there are two possibilities. For small values of B > 0, both components 
continue to be monotone increasing in r. But (for instance) when A + < B < A_, the 
component f+(r) has /+(0) = but it approaches its limiting value from above, and hence 
cannot be either increasing or decreasing in r. This is a phenomenon which is not familiar 
in the GL theory, although non-monotone vortex profiles have been found in the context of 
confined BEC (see jKoP] ). 



Parts (i) and (ii) of Theorem 1.2 are proven using the second variation of energy around 



an equivariant solution. Part (iii) is obtained via a compactness argument based on the 



uniform asymptotic estimate from Theorem 3.2 



The questions addressed in this paper are motivated by analogous results for the classical 
(scalar complex) Ginzburg-Landau equations, 

- Am + (H 2 - l)u = 0. (11) 



The existence, uniqueness and strict monotonicity of equivariant solutions to (11) were 
proven in [HHJ, and the asymptotic expansion of the radial solutions as r — > oo was de- 
rived in |GEQ . Results for a system of two Ginzburg-Landau equations ([I]) have been 



recently been obtained in [AB], [AB Mlj . |ABM2] in the special case of "balanced" coeffi- 
cients, A + = A^ and t + = £_ = l/y/2. For the balanced case, existence, uniqueness and 
monotonicity of radial solutions are proven in [ABMlj, and questions of minimality, stability 
and bifurcation are treated in [ABM2J . While several of the methods employed in these pa- 
pers may be generalized to the general case (with hypotheses (H)), the balanced coefficient 
case provides several simplifications which are not valid in the general case. For instance, in 
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the balanced case with n + = n_, the unique equivariant solution to ([!]) takes the simple form 
$ = f^=, ^=J where it solves the classical (complex scalar) Ginzburg-Landau equation (11), 
and thus the full treatment of the system of equations ([!]) may be reduced to dealing with 
a single equation. As discussed above, the general case provides more variety of possible 
behaviors for symmetric vortex solutions, and requires enhanced methods which apply to 
systems of equations. 

A forthcoming paper will consider the stability and minimality of the symmetric vortex 
solutions. 

Two- component Bose-Einstein condensates 

In a model of a two-compenent BEC from Eto et al. [EK NTTj . we consider a pair of complex 
wave functions $ = (tpi, ip 2 ) £ C 2 , defined in the sample domain Q C M 2 . The energy of the 
configuration is defined as 



h 2 ,_ ,„ h 2 ,_ ,0 1 



|Vy?i| 2 + 7, — |Vy?2| 2 + - (fi-i^il 4 + 92^2^ + 2g 12 \ip 1 \ 2 \ip 2 \ 2 ) 



2mi 2m 2 2 



dx. 



where mx,m 2 > are the masses, and the coupling constants satisfy the positivity condition 
9\92 ~ 9\ 2 > 0. The Gross-Pitaevskii equations govern the dynamics of the condensate, 



h 2 

ihd t <fi = -- — A<p x + 5-i|y9i|Vi + g 12 \<fi 2 \ 2 (Pi, 
2m 1 

h 2 

ihd t Lp 2 = -- — Ay? 2 + gi 2 \(pi\ 2 ¥>2 + g 2 \<P2\ 2 <P2- 

2m? 



(12) 



A stationary equation of the desired form is obtained by considering standing wave solutions, 
(fi(x,t) = e~ ltHt / h Ui(x), i — 1,2, where /Zj represent the chemical potentials: 

h 2 



-Acpx + gx\<pi\ 2 (pi + g 12 \(f 2 \ 2 (pi = fii<p 



2mi 

h 2 

-A(p 2 + gi 2 \^i\ 2 f 2 + 92\^2\ 2 ^2 = H2V2 



2m 2 

In the variational formulation of the stationary problem, the chemical potentials represent 
Lagrange multipliers, which arise because of the constraints on the masses of the two con- 
densate species, 

\ipi\ 2 dx = / \ui\ 2 dx = Ni, i = 1,2. 
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By a rescaling of the dependent variables, ip + = ^^u\, ■0_ = \J^~ 2 u ii we may eliminate 
the masses from the equations, and we obtain the system 

- e 2 A^ + + L4 + (|^ + | 2 - t 2 + ) + B{\^\ 2 - t 2 _ = 0, 
-e 2 A^_ + [A_(|^_| 2 - t 2 ) + 5(|^+| 2 - t 2 ,)]^ = 0, 

with e 2 = -^=, = 2*0!, A_ = ^ 2 5 = 12 , and 

2 _ A*lff2 ~ A*2ffl2 / mj 2 _ A*201 - ^1012 



0102 - #12 V m l ' #102 - 5 , ?2 




As is typical for GL equations, by blowing up around the core of a vortex at scale e, we 
obtain an entire solution of (jl]) in all of 1R 2 . 

In a more physically appropriate model for a two-component BEC |KTUj . the Laplacian 



in the Gross-Pitaevskii system (12) should be replaced by the Hamiltonians 



2m; 

with harmonic trapping potentials V^ trap = c 2 |x| 2 /2, c; constant i = 1,2; angular momentum 
operator L z ; and (constant) angular speed Q. While this is an essential step, both in modeling 
the confinement of the condensate and in describing the onset of vortices in the sample, these 
terms will not affect the general form of the blow-up equations which describe the vortex 
profiles at length scale e in the condensate. Indeed, the momentum operator plays much the 
same role as the magnetic vector potential in the GL model of superconductivity, and for 
rotations which are of moderate strength in e, u <C e -1 , the analysis of Proposition 3.12 of 
[SSJ may be used to derive ([I]) in limit e — > after rescaling. 

2 Existence and uniqueness 

We begin with an a priori bound on solutions of the system (jl) in bounded domains. Let 

a _£? 

A, > denote the smallest eigenvalue of the matrix J" . 

to I B A_ 

Proposition 2.1. Let Q C M 2 be a bounded smooth domain, and ^ = (ip+, ipJ) be a solution 
of (jlj in Vt, satisfying 

\^(x)\ 2 = \tp+(x)\ 2 + \ip-(x)\ 2 < t\ +tl for all xe dVl. (13) 

Then \^(x)\ 2 < minj 2 ^, t\ + t 2 _ } for all x E Q, where M = m&x{A + t 2 + + Bt 2 _ , A_t 2 _ + Bt\}. 
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We note that the bound on H^Hoo depends only on the coefficients in the energy, and in 
particular it is independent of the domain Q. Furthermore, the quantity M > whenever 
B 2 < A + A_. 

Proof. Let V(x) = \^\ 2 . We then calculate: 

l -AV = W+) + W-) + |Vvl/| 2 

> (^+,W>+) + (</>-, W>_) 

= [A + {\^ + \ 2 - t\) + 5(|^| 2 - t 2 _)]|^+| 2 + H_(|^_| a - t 2 _) + S(|^ + | 2 - i 2 )]|V-| 2 
= [A + |V + | 4 + 2 J E# + | 2 |^_| 2 + A_|V_| 4 - (A+t 2 , + fit 2 )|^ + | 2 - (A_t 2 + 5t 2 + )|^_| 2 ] 

> [A s (|^+| 4 + |^_| 4 ) - (A + t 2 + + fit 2 )|V + | 2 - (A_t 2 + BtlM_\ 2 } 

> ( ^|^|4_ M |^| 



Now let f = V— ^r~, which is then a subsolution of the following boundary value problem: 

-Av + a(x)v < in fi, 
< + t 2 - 



an A s 
with a(x) = \ S V > 0. 

If t\ + t 2 — 2 ¥- < 0, by maximum principle t> < on Q, i.e. V < 2 r^-. On the other 

i As 

hand, if t 2 , + t 2 _ — ^ > 0, by maximum principle again rnaxf can be attained on <9f2, i.e. 
max?; < t 2 + t 2 _ — and v < t 2 + t 2 _ — , i.e. V < t 2 + 1 2 _. We thus conclude that 



A s A s 



V <min|— ,i+ + t_ >, i.e. |^| 2 < min j — ,t 2 + t 2 J>. □ 
We may now prove the existence and uniqueness theorem. 



Proof of Proposition 1.1 Without loss of generality, we assume n± > 0, as negative degrees 
may be obtained by complex conjugation of ip±. To obtain the existence we consider the 
problem defined in the ball Br, R > 0, 

-fi ~ ~f± + %f± + [A ± {fl-t 2 ± ) + B{fl-tl)] f ± = 0, for < r < R, 
/±(0) = if n ± ^0; /i(0) = if n± = 0. 
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The existence of such a solution follows easily by minimization of the energy 



If £ 



,i=± 



rit 1 



(fir + i±f? + =;iMf? - *?) a + 25(4 - - a 



r 2 '* 2' 

(16) 

over Sobolev functions satisfying the appropriate boundary conditions at r = and r = R. 



Denote [fR,+ {r), fii,-{r)) as any solution of ([THJ . Since E^ +> n _(|/+|, |/_|) = E* +> „_(/+, /_), 



and there exists a minimizer which is nonnegative, f± > in [0, ii]. 



By Proposition 2.1 , the solution pair is bounded, + /#_(?") < A 2 := minj 2 ^, i 2 + 

£ 2 }, for any solution to (|3]), uniformly in i?. By L p estimates (applied to each equation 
separately), we may conclude that f^± is uniformly bounded in Wf£ for any p, and by 
elliptic regularity the family is bounded in Cf oc for any k. Thus, there exists a subsequence 
R n -> oo for which the solution [fn t+ {r), f R ,-{r)) [/«>,+ (>"), foo,-{r)] in Cf oc ([0, oo)), and 
the limit functions [foo,+{r), foo-(r)} give (weak) solutions to the ODE on (0, oo) with the 
same boundary condition at r = and satisfying the uniform bound + (r) + _(r) < A 
for all r G [0, oo). By uniform convergence, we also have /oo,±( r ) > for all r G [0, oo). 

To determine the behavior near r = 0, we write the ODE in the form: 

-£-^ + =|/ ± = fc(r), 

with g± smooth, and <7±(0) = for n± ^ 0. By the Frobenius theory (see Appendix 3 of [BCJ 
or Lemma 5.9 of |ABM2j ). we deduce that the behavior foo,± ~ rU± near r = 0, which is (J6|. 
As a consequence (see Theorem 3.1 of |BCj ). we may conclude that if)±(x) = foo,±{ r ) e% is 
regular at r = and solves (nj) in M 2 . 

Next we must verify that /oo,± attains the desired asymptotic conditions as r — > oo. We 
first claim that 

/•oo 

/ [(^, + (r)-4) 2 + (/^(r)-t 2 ) 2 ]rdr<oo. (17) 
Jo 

For this purpose, we derive a Pohozaev identity: we multiply the equation of fn t ± by r 2 f' R ± (r) 
and integrate by parts with respect to r G (0, R). Adding the two resulting identities gives: 

[Rfk + (R)? + [Rfk-( R )? 

+ T {Mfi+ - 4) 2 + Mfl- - V + 2B (fl + - tWl- - = 
Jo 

nitl + nV, (18) 
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which is the Pohozaev identity. By the uniform convergence f^± — >■ foo,± on [0, Rq] for any 
Rq > we have 



A, / [(A + (r)-^) 2 + (/^,_(r)-e) 2 ]rrfr 

< / {A, - t\f + AMI,- - tiT + 2B{£ )+ - tl)(fl_ - £)} rdr 
./o 



< niti + n 2 _t 2 _ 



for any fixed _R > 0. Letting R — > oo we establish the condition (17). 



Next, we prove (|7|. In the remainder of the proof, we write f± = foo,±- Let ^i(r) be 
a smooth function with ^i(r) = 1 for < r < 1 and T)x(r) = when r > 2, and set 
? 7i?( r ') := Vi( r /R)- We then multiply the equation for / ± by r(t± — f±)i]R an d integrate by 
parts to obtain: 



f f±(t± - f±)vR + Mfi - tl)f±(t± - /±) + B{f 2 T - t%)f ± {t ± - u) 



/OO 1 POO 

(f±) 2 VRr dr + - J ((t ± -f ± ) 2 )'r]' R rdr 

POO /* oc 

/ (/±) 2 VR r dr - - (t± - f±f A r T] R rdr, 
Jo 1 Jo 



r dr 



where A r f := -(r/'(r))' is the Laplacian for radial functions. We now estimate each term, 
using (fX7b. We note that, as /± > 0, we have {t\- flf = (t± + f±) 2 {t±- f±) 2 >t 2 ± {t±-f±) 2 , 



and hence \\t± — /±||l 2 < i il^± — /±iU 2 < oo by (17). Thus, we obtain the bounds (uniform 
in R), 

r < poo 

(t ± -f ± ) 2 \A rVR \rdr<— / ( t± -f ± ) 2 rdr^0, 
o K Jo 



oo 2 

o r2 



< A 



ll*±-/±IU" <Ci, 



2 .2 



< A ± A||^-4|| L2 ||/ ± -t ± || L2 <c 2 



-f±(t± - f±)r]Rrdr 

/■oo 

/ A ± {fl-t 2 ± )f ± {t ± -f ± )rdr 
Jo 

and similarly for the B term in the integral. We thus conclude that 

POO 

/ (f' ± ) 2 VRrdr<C 3 
Jo 

is bounded independently of R, and the finiteness of the integral ^ follows from Fatou's 
lemma. 
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The limit statement /±(r) — > t± as r — > oo now follows from Proposition 2.1 and the 
bounds (17) and 0. Indeed, for any < rj < r 2 , 



(/±(r- 2 ) 



< 2- 



A 

i 

Vi 



''2 



(/ ± (r)r + (/l-4) 2 



oo. 



and hence (f±(r) —t±) 2 is Cauchy as r 

For the existence part of the theorem, it remains to verify ([8]), we note that by ^ there 
must exist a sequence R n — > oo such that R n f'±(R n ) —¥ 0. Applying (18) with i? = R n , and 
passing to the n — > oo limit gives ([8]). This completes the existence part of Proposition |1.1| 

Uniqueness may be proven in (at least) two different ways. The first employs the method 
of Brezis and Oswald |BOj . which we now sketch. Let [n+, n_] G Z 2 to be given, and suppose 
[/+, /_] and [g + , gJ\ are two solutions of (|3]). Dividing the equations for f± by f± (and 
similarly for g±), multiplying by f\ — g\ and integrating by parts, we obtain: 



and similarly, 



[A + (fl-9l) 2 + B(fl-9l)(fl-gi)] rdr 

(/+ - 9 2 +)rdr 



A r (U) A r (g, 



0H 



(/;) 2 -2^/x+fw 2 



+ 



(9' + ) 2 -2 g -±f + g' + + 9 -±(f + r 
J+ J+ 



f i f+ i 

U - —9+ 
9+ 



Q > _9±f 

y+ f j+ 



rdr 
rdr 

rdr, 



f- - —sL 

9- 



t. - jL 



rdr 



[A-(f--9l) 2 + B(fl-gl)(f_-gl)} rdr. 



Adding the above identities together we obtain that 



< 



/; - -9\ 

g+ 



g+ £ , 



9 + ~ T-f + 



f- ~ -9- 

g- 



rdr 



POO 

/ {A + (fl ~ 9lf + AM 2 ~ 9 2 -) 2 + mil - 9l)U 2 - 9 2 -)}rdr < 0, 
Jo 



(19) 
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as the quadratic form on the right-hand side is positive definite. Therefore, f± — g± = 0, i.e. 
f± = g± for all r e (0, oo), and this completes the first proof of uniqueness. 

For a second proof, we note that the system of equations for f± may be written in the 
form 

-AU + d u G(r,f 2 + J 2 _)U = 0A 
-A/_ + d v G(r, ft, /?)/_ = 0, 1 
f±(r)>0, J 

with 

n 2 ?i 2 1 

G(r , u, v) = ^u+^v+- (A + (u - t\) 2 + 2B(u - £)(t; - t 2 _) + A_(v - t 2 _) 2 ) . 

As G is convex in (u,v), we may apply Theorem 4.1 of |ABM1] (together with the limiting 
argument given in that paper, to apply the theorem to the semi-infinite interval) to conclude 
uniqueness. □ 



3 Asymptotics 

We derive the asymptotic behavior of solutions f±(r) as r — > oo by means of the sub- 
and super-solutions method. The basic tool we require is the following useful and original 
comparison lemma: 

Lemma 3.1. Let A, B,C,D be bounded functions on [R, oo). 

(A) Assume A, D > 0, B, C < and 4AD - (1 + C) 2 > in [R, oo). Then, if u, v are 
radial solutions of the following problem 

( -Au + Au + Mv <0 in [R,oo), 
\ -Av + Cu + Bv < in [R,oo), 

withu(R) < 0,v(R) <0,u,v are bounded in [R,oo) with {u') 2 rdr < oo and (v') 2 rdr < 
oo, we have that u < and v < in [R, oo). 

(B) Assume A, D > 0, B, C > and 4AD — (B + C) 2 > in [R, oo). Then, if u, v are 
radial solutions of the following problem 

( -Au + Au + Mv < in [R,oo), 
\-Av + Cu + Bv >0 in [R,oo), 

with u(R) < < v(R), u,v are bounded in [R, oo) with \u') 2 rdr < oo and f R {v') 2 rdr < 
oo, we have that u < < v in [R, oo). 
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Proof. To verify (A), multiply the respective equations by u + = max(u, 0) and v + = 
max(f,0) and integrate by parts in [R, T] for VT > R > 0. Since u and v are radial 
solutions, we have 

Au ■ u + dr = — I (ru')'u + dr = — ^u'u + r^j + I [(u + )'] 2 rdr. 



hence, 



i.e. 



T pT pT 

Au ■ u + dr = — {ru')'u + dr = — ( uu + r\ R \ + / 

R JR ^ ' J R 

u'u + r\ J + / [(u + )'frdr + / \k{u + f + B(v + + v-)u + ]rdr < 0, 

J Fi J Ft 



\{u + )'} 2 rdr+ I [K(u + f + Mu + v + + Mu + v~]rdr 

R JR 

< u'(T)u + (T)T - u(R)u + (R)R 

< u'{T)u + {T)T , (20) 

since u\R)u + {R)R = by u(R) < 0,u + {R) = 0. 
Now we make a claim: 

Claim A: There exists T n — > oo such that u'(T n )T n — > 0. 

Argument by contradiction. If not, there exist co,Tq > with u'(r)r > cq for all r > To. 

2 

It follows that u'(r) > a , then we obtain that [u'(r)] 2 r > ^ ^ L 1 , which is a contradiction. 

Combine boundness of u in [R, oo) and the result of above claim, we take T — > oo on the 
both sides of (20), we get that 

) poo 

[(u + )'] 2 rdr + / \A(u + ) 2 + Mu + v + + Mu + v~]rdr < . (21) 

R J R 

Similarly, we have the inequality for v + : 

[(v + )'] 2 rdr + / p(v + ) 2 + Cu + v + + £u~v + ]rdr < . (22) 

R J R 

Since Mu + v~ > 0, Cu~v + > with u + = max(-u,0) > 0, v + = max(f,0) > 0, v~ = 
min(f , 0) < 0, we get that 

POO 

\(u + )'} 2 rdr+ j [A(u + ) 2 + Mu + v + }rdr 

/>' J R 

POO POO 

< / [(u + )'] 2 rdr+ / [A(u + ) 2 + Mu + v + + Mu + v~]rdr 

J R J R 

<o, 
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and 



[(v + )'] 2 rdr + / [B(v + ) 2 + Cu + v + }rdr 

11 J R 

OO POO 



< / \{v + )'] 2 rdr + 

J R J H 

< 0. 



[V 



+ Cu + v + + Cm v + }rdr 



Therefore, we deduce that 



R 



{[(« 



+ V12 



\/l2 



A(w H 



) 2 }rrfr < . 



(23) 



Note that the matrix associated to the quadratic form is positive definite in [R, oo) by 
hypothesis, so there exists a function A + > with 



A(u + ) 2 + (B + C)u + v + + B(v + ) 2 > A+ [(u + ) 2 + (v + ) 2 ] > 



In consequence, 



< 



{[(« 



+\/l2 



+ [(« 



+ \/l2 



A + [(w + ) 2 + (t; + ) 2 ]}rrfr<0 



which implies that u + = = u + in [i?, oo). Therefore u(r) < and u(r) < in [R, oo). 

To prove (B) we multiply the first inequality by u + , but multiply the second inequality 
by v~ = min(t>, 0) < 0. Then, integrate by parts, we get that 

T pT 

[(u + )'] 2 rdr + / [A(u + ) 2 + Mu + v + + Mu + v~}rdr < u'(T)u + (T)T 

R J R 

and 

/ \(v~)'] 2 rdr + / [Cm + w~ + Cu~v~ + D(w") 2 ]rrfr < v'(T)v~(T)T 

J R J R 

by f(i?) > and v~(R) = 0. By the same argument as in Claim A above, there exists 
T n — > oo such that v'(T n )T n — > 0. Hence, passing to the limit T n — » oo, we have 



[(w^)'] 2 ^ + / [Cm + w~ + Cw~*T + D(w") 2 ]rrfr < 0. (24) 

' R J R 

Since Mu + v + > 0, Cm^w - > with u + = max(w, 0) > 0, u~ = min(w, 0) < 0, v + = 
max(f,0) > 0, v ~ — win(v, 0) < 0, we can drop Mu + v + and Cu~v~ two terms without 



affecting the inequalities (21) and (24). We add (21) and (24) together, and note that the 



matrix associated to the quadratic form is positive definite in [R, oo) by hypothesis, so there 
exists a function A + > with 



u +f + (B + C) u + V - + B(v~y > \ + [(u + Y + [v-f] > . 
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R 



R 



' } rdr < . 



Consequently, 



0< / {[{u + )f + [{v-)f + \ + [(u + ) 2 + (v-) 2 }}rdr <0 , (25) 

lR 



which implies that u + = = v ~ in [R, oo). Therefore u(r) < and v (r) > in [R, oo). □ 

After this preliminary, we introduce the asymptotics of the radial solution at oo. 

Theorem 3.2. Let [/+,/-] be the solution of ^ with degree pair [n + ,n_] at oo, then we 
have 

f±=t± + ^ + ^ + 0(r- 6 ) as r^oo, 



with 



and 



a± 



1 Bn\ - A T n 2 ± 
2{A + A_ -B 2 )t± 



(26) 



_A 2 f (8n 2 ± + n%)t 2 T - BA T (2n 2 ± + 8)n|t| - &BA±n%t 2 ± + B 2 (8t 2 ± n 2 ± + n%t%)t 2 ± 
±_ 8(A+A_ - S 2 ) 2 4 1| • ( 27 ) 

More specifically, let A± > and i?o > so that A + A^ — B 2 > 0. Then, there exist positive 
constants C\,C2,R such that 



f±(r) 



a± 



< 



Ci 



2a 



± 



< 



C 9 



(28) 
(29) 



hold for all r > R and all B, \B\ < Bo- 
Proof. We treat the cases B > and B < separately. We start with B > 0. 
Step 1: Construction of a sub-supsolution pair. Let 



£+ + a H 
t_ + a_ 



R 2 


+ b+ 


i? 4 


i? 6 








R 2 


+ 6_ 


i? 4 


R 6 




r 4 





(30) 
(31) 
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where a±, b±,c + ,c_ and R are to be chosen so that 

w' 1 



LHS(w+) 
LHS(w_) 



-w" 



-w 



u/_ 1 
h 



r 



[A + {w 2 + - t\) + B(w 2 _ - t 2 _)]w + > 0, 
[A-(wi - t 2 _) + - tV)]w_ < 0, 



(32) 
(33) 



for all r > R, and < w+(i?), /_(#) > w_(i?). 



Expanding (32) and (33), we obtain terms which are polynomials in even powers of r x , 

(34) 



LHS(^ + ) = ]Tm+ - , 
fc=i ^ ' 

9 /f>\ 2k 

LHS(^ = £M 2 - - , 



(35) 



where M* = M±(A±, B, R, 

o±> b±, c+, c_) is a polynomial in its arguments. The coefficients 
are quite intricate, but may be explicitly calculated with the aid of Maple software. 

First, we choose a± in order to force the lowest order coefficients M 2 to vanish: that is, 

M£ = [1 + 2{A±t±a± + Bt T a T )R 2 ]t± = 



is achieved by choosing a± as in (26). Similarly, we fix the values of the coefficients b± in 
order that the next term 



Mi 



{-3a±R 2 + 2R A (A±b±t 2 ± + Bb T t±t T + Ba±a T t T ) + R 4 t±(3A±a 2 ± + B 2 )]t ± = 



which determines the values for b± given in (27). 

Again using Maple, the values of a±,b± may then be substituted into the expansions of 
(32) and (33), and the expressions for M 2k may be viewed as functions of R. The exact form 
of the coefficients M 2k is very complex, but they are all rational expressions of the coefficients 
A±, B, t±, n±, and as we will choose R large, we are only interested in the leading order of 
each. We obtain: 

M+ = 2(Bc_t_ + A + c + t + )t + + 0(R~ 6 ), 



Ml 



2(Bc + t + + A_c_t_)t_ + 0(R~ 6 ), 



± 



0{R~ 



M 



10 



0(R~ 



M 
M 



12 



12 



(A + c 2 + + Bc 2 _)t + + 2(Bc^ + A + c + t+)c+ + 0(R~ 
(A_c 2 _ + Bc\)t- + 2(Bc + t + + A_c_t-)c_ + 0{K 



M 



ii 



= 0(iT 2 ), M± = 0(iT 4 ) 
M+ = (A + c 2 + + Bc 2 _)c + , 



M 



18 



(A_cl + Bc 2 + )c_, 
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here 0(R~ n ) denotes terms which are small as R — > oo uniformly for \B\ < B . 

As is the leading order term, we choose c+, c_ in order that it gives the correct signs, 
and so that it dominates the other terms for r > R. To this end, we let [c+, c_] be the unique 
solutions to the following system: 

A + t + c + + Bt_c_ = 1 , 
Bt+c+ + A_i_c_ = -1 . 



We note that, 

c+ 

Let c + = 5c + , c_ 



A_+B 



>0, 



A, + 5 



< . 



- B 2 )t + ' (A+A_ - B 2 )t. 

5c_ with < 5 < 1 to be chosen later, hence c + > 0, c_ < 0, and 



M+ = 25t + > 0, 



-25t_ < 0. 



By choosing R sufficiently large, we obtain that 



I Mo 



IM^I, 



and so we have that 



Mf 



10 



14 



20' 



16 



\M, 



6 h 



Also, 



(36) 



(37) 



(38) 



for Vr > i2. (39) 



M+| = \5\A + c% + Bci)f + + 25 z c + 
= 5 2 \{A + c 2 + + Bc 2 _)t + + 2c + \ 
< Cl (A±,B,t±)5 2 



5 K|, 



for R sufficiently large, provided we choose 5 < 
Meanwhile, 

= |5 2 (A+5^ + B5 2 _ )t+ + 25 2 g+| 



M- 



121 



= 5 2 \(A_c 2 _ +25_ 
<c 2 (A ±)J B,t ± ),5 2 



< -<ft_ = - 
~5 5 



M, 



6 h 



(40) 



(41) 
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for R sufficiently large, provided we choose 5 < Therefore, if we choose 5 < mm{^- } %r} 



5c2 

it yields that |M^| < ||M^|. So we can deduce that 



M- 



R 



12 



< 



mi 



R 



(42) 



Finally, we note that Mf^ has the same sign as Mg , and so it contributes with the desired 
sign. 



By (|34|)-(|35j) and (|37j), (|38j), (J39J), (|42j), we deduce that 

9 / n\ 2k -, 



^=1 



> --M, 
5 



6 



-M? — 
5 6 V r 



6 / -R 
5 + \ r 



6 



6 



>-5t+f-l >U. 



(43) 



fc=i 



2A- 



< 



I AG 



4m 6 



7? 

r 

/? 

r 



+ -|M r : 



< 



:St. 



6 



< 0, 



5t_ 



R 

r 

R 

r 
R 

r 



+ m: 



+ 5 M e 



for all r > R. Thus, iu + , u>_ satisfy (32) and (33), as desired. 



R 

r 

R 

r 



(44) 



Now we consider w+ and w_ at r = R. Combine (30), (31) and (26) we have 



t + + c + -0{R 



-2\ 



W 



t- + c_ —0(R 2 ). Since f±(r) — > t± as r — > oo and c + > 0,c_ < 0, 



we deduce that w + (R) > f+{R) and w_(R) < f-(R) for R > Ro sufficiently large. This 
completes Step 1. 



Step 2: We will apply Lemma 3.1 to show that w+(r) is a super-solution to / + -equation 
of (|3| and w_(r) is a sub-solution to /--equation of Let u = f + — w + , v — /_ — 
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together with equations d3|), (32) and (33), we have that 



and 



—A r u = —A r f + + A r w^ 



(wl 



A + (t 2 + -f 2 + )f + + B(t 2 _-f 2 )f + + A, 
A + (w% - fl) - A + t 2 Jw + - f + ) + Bt 2 _ (f, 



r + )w+ + B(wi - ti)w+ 

w + ) + B(wlw + -flf + ) 



[A + t 2 + + Bti - A + (w% + w + f + + fl)\u + B(wiw + - fif + ) 
[A + t\ + B(t 2 _ - f 2 ) - A + {w 2 + + w + f + + f 2 + )]u - B(w_ + f_)w + v 



-A r V 



-A r /_ + A r w_ 

A-(tl - /!)/_ + B(t 2 + - fl)f_ + A_{w 2 _ - t 2 _)w_ + B{w\ - f*>_ 
A-(wl - /!) - A_t 2 («;_ - /_) + - w_) + S(tZ)5.ul- - /+/ 

[y4 _t 2 + £?(^ - /2) _ A_(^ 2 + W_f_ + / 2 )]^ - + f + )w_U , 



and thus (32) and (33) imply that 

- A r u + [A + (w 2 + + w+f+ + / 2 ) - + B(f 2 - t 2 _)}u + B(w_ + f-)w+v < 

- A r v + 5(w+ + f + )w_u + [A_(w 2 _ + + f 2 ) - A_t 2 _ + B(fl - t\)\v > 
As in Lemma 3.1| we have a system on [R, oo) of the form: 

-A r u + A(r)u + M(r)v < 0, u(i2) < 0, 
- A r v + C(r> + B{r)v > 0, > 0, 

with 

A(r) = A + (w 2 + + w + f + + f 2 + ) - A + t\ + 5(/ 2 - t 2 _ ), 
B(r) = + f-)w+, C(r) = + f+)w_, 

D(r) = A_(«; 2 + w_f- + /!) - A_t 2 _ + 5(/ 2 - £). 



(45) 
(46) 



Now, we check the hypothesis of Lemma 3.1 with the uniform convergence of f±(r) as r — > oo: 

A(r) ->■ B(r) -> 2£t+t_, 

C(r) -> 2Bt + t_, D(r) ->■ 2A_t 2 , 

and 



4AE 



+ C) 2 -> 4 • 2A + t; • 2A_r - (4B*+t-) = M>t + t_(A + A- - B 2 ) > 



for i? sufficiently large. All conditions of Lemma 3.1 are satisfied, it yields that u(r) < 
0, v(r) > in [R, oo), i.e. f+(r) < w+(r), f-(r) > w_(r) in [R, oo). 
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Step 4: We repeat the above procedures with the roles ± of w + and w_ under the 
condition B > 0. We may obtain c + , c_, R so that /+(r) > w + (r), < w_(r) in [R, oo). 

Combine the results from above steps, we get that 



w+( r ) < /+( r ) < w +( r )? w_(r) < /_(r) < w-(r), for Vr > R, 



(47) 



in the case B > 0. 

Step 5: Now we consider the case 5 < 0. This is similar to the previous cases, but we 



must set up part (A) of the Comparison Lemma 3.1 Firstly we construct a super-solution 
pair, with w + as in (30), u)_ as following, 

t_ + 



R 2 



R 4 



R 6 



with a±,b±, c± and R are to be chosen so that 



LHS(u>+; 
LHS(w_; 



-// 
-w. 



-w" 



w' 



+ —w+ + [A + (w 2 + - t%) + B(wi - t 2 _)}w + > 
+ + [A_(w 2 _ - t 2 _) + B(w\ - t 2 + )]w_ > 



(48) 

(49) 
(50) 



for all r>R, and f+(R) < w+{R), f-{R) < w_{R). 



Using Maple software, we expand (49) and (50), which is a polynomial in even power of 
r -1 , as in (34) and 

9 



LHS(^ = £M- fe (^) 
k=i ^ ' 



2/v 



(51) 



respectively, where M^ k = M£ e (A±, B, R, a±, 6±, c±) is a polynomial in its arguments. Sim- 
ilarly, we choose a±,&± as in (26) and (27) separately so that may be zero, arriving 
at the same formulas in Step 1. Then, repeat the same procedure as in Step 1, we get the 
exact forms for M 6 , M 12 an d M 18 as the same as in Step 1 but with c_ replaced by c_, and 
M^ 1 , and as the same as in Step 1. 

As is still the leading order term, we choose c± in order that it gives the correct sign, 
and so that it dominates the other terms for all r > R. To do this, we take [c + ,c_] to be 
the unique solutions of the linear system: 

A + t + c+ + Bt-d- = 1, 



Bt + c+ + A_t_c_ = 1, 



that is: 



A_ -B 



(A + A_ - B*)t+ 



> 0, 



A.—B 



(A + A_ - B 2 )t. 



> 0. 
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Let c± = Sc± with < 5 < 1 to be chosen later, hence c± > and 



Mt = 26t± > 0. 



(52) 



By choosing R large enough, we can also get (38), (39) and 



\ Ad- 



vil 



= \S 2 (A + c 2 + + Bc 2 _)t + + 25 2 
= 5 2 \(A + c 2 + + Bc 2 _)t + + 2c^ 
<c 5 (A ± ,B,t ± )5 2 



<\m« 



-Ma 



(53) 



it j 



which implies that 5 < 5 
Meanwhile, 



|M- 



121 



= \5 2 (A_c 2 _ + Bc 2 + )t_ + 2S 2 c. 
= 5 2 \(A_c 2 _+Bc 2 + )t^+2c_\ 
<c 6 (A±,B,t±)5 2 



which implies that 5 < 



2t- 
5c 6 



(54) 

Therefore, if we choose 5 < min{|^,|^-}, it follows that 



\Mf 2 \ < \\M±\. And for Mf % , we have 



\M±\ = 5 3 (A ± c 2 ± + Bcl)c ± = 0(S 3 ) 



uniformly in A±, B, and thus it can still be controlled by M 6 . In particular, taking S smaller 

ii 

5 I 



if necessary, we have | M^ 8 | < 1 1 M§ \ 



By (|34j), (|5ip and (|52j), (|38j) - (|39j) and (|42j), we deduce that 

9 



J2 M 

k=i 



2k 



~ 5 6 1 r 



1 I R 

-Mr 



'6 



> 0, 



(55) 



for all r > R. 



Now we consider w± at r = R. Combine (30), (48) and (26), we have w± = t± + c± 



0(R 2 ). Since f± — > t± as r — > oo and c± > 0, we obtain that w±(R) > f±(R) for R > Rq 
sufficiently large. This finishes Step 5. 
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Step 6: We will apply Lemma 3.1 to show that w±{r) is supsolution to (fl]). Let 

/+ — vo + , v — /_ — w_, we have that 

—A r u = —A r f + + A r w + 

= L4 + (4 - fl) + 5(t 2 - / 2 )]/ + + [A + {wl - t\) + B{vo 2 _ - ti))w + 
= -A+UI ~ <) + ^V+(/+ " tD + ) + - tD + ) + B(w 2 _ w + - f_f + ) 

= [-A + (fl + f + w + + w 2 + ) + A A + Bt - ~ B f-\ u ~ B (™~ + f-)™+ v 
= [A + t 2 + + B(£ - f 2 ) - A + (fl + f + w + + w 2 + )]u - B(w_ + f-)w + v, 



u 



and 



-A r v 



-A r /_ + A r w_ 

[A^tl - f 2 ) + B(tl - fl)]f- + [A-(t/£ - £) + B(tD* - £)]«;_ 
-A_(/! - u£) + A_£(/_ - u)_) + - tD_) + 5(«; 2 «;_ - / 2 /_ 

[-A_(/ 2 + + w 2 ) + A_t 2 _ + St 2 - 5/ 2 ]i; - B{w + + 

-S(t5 + + f + )w-u + + + ^ 2 ) + A_t 2 + B(t 2 + - f 2 + )]v, 



thus by (|49j)-(|50j) we have that 

-A r u + [A + {f + + + w 2 + ) + - - A + t2 +] u + + f-)^+v < 0, 

-A r i; + B(vo + + /+)«>-« + L4_(/! + + w 2 ) + B(fl - t 2 + ) - A^]v < 0, 

and u(R) = f+(R) — w + (R) < 0, v(R) = f-(R) — w~(R) < 0. Therefore, we have the exact 



(56) 

(57) 
(58) 



forms of the inequality system as in the part (A) of Lemma 3.1 with 

A(r) = A + (f 2 + + f + w + + wl) + B(f 2 -tl)- A + t 2 + , 

B(r) = -B(w- + C(r) = -B(w+ + /+)«>_, 

D(r) = A_(/ 2 + + «; 2 ) + B(/ 2 - 4) - A_i 2 _. 

With the uniform convergence of /±(r) as r — >■ oo, we obtain that 

A(r) -> 2A + t 2 + , B(r) -> -2Bt + t_, 

C(r) -> -2Bt+t-, D(r) -> 2A_t 2 

and 

4AD - (1 + C) 2 ->■ 4 • 2A + t^ • 2A_t 2 - (4BM-) 2 

= 16i+t 2 (A + A_ - fi 2 ) > for i? large enough. 
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Hence we can simply apply part (A) of Lemma 3.1 to get that u(r) < 0, v(r) < in [R, oo), 
i- e - f±( r ) < w±(r) in [R, oo). 

Step 7: We repeat the same procedure as in Step 5 and Step 6, but with a pair of 
sub-solutions w± to ([T| instead of sup-solutions w±. Let w + be defined as in the following 

p2 p4 p6 

it , it rt . . 

w + = t + + a + — + b + — + c + — , (59) 



and w_ defined as in (31 ), with a±, b±,c± and R are to be chosen so that 

LHS(w ± ) = -w!i - ^ + ^w ± + [A±{wl - t 2 ± ) + B{wl - t%)]w ± < 0, (60) 

-Sc± with the same S as in Step 



for all r > R, and w±(R) < f±(R). And we choose c ± 
5 so that M§ is still the leading terms in (34) and (51) with which w± replaced by w ± . 
Therefore, we can simply get that 



LHS(w ± ) < Ml 



R 



r 



+ j:\Mt\ 



Mi - - l -M± 



:M± 



<> 



--5t (- 

5 ± \ T 



< 0. 



(61) 



5 ~ \r 

Next, with u = /+ — w + , v = /_ — w_, we repeat the similar process as in Step 6 and 



apply part (A) of Lemma 3.1 again, we can obtain that w±(r) < f±(r) in [R, oo). 



Step 8: Combing all the results from above steps, we have that 

w±( r ) < f±( r ) < w±(r) in [R, oo). 

for both cases of B. 

From the sub- and super-solution argument, we have that 



(62) 



< 



with a± = a±R 2 , b± = b±R 4 , a± and b± defined as same as in (26) and (27) respectively, and 
C\ = mm{c±R 6 } with c+ = min{|c+|, |c + |}, c_ = min{|c_|, |c_|}. 



Next we want to show that \f'±{r) + < ^ for some constant C^- Following the 

: f± ~ w ± + % (f° r convenience, we drop 



idea in [CEQ] , let W ± (r) = / d 
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bar and underline of w± and c± as in the formulas of w± shown in previous steps), hence 
W±{r) = 0(r -6 ). Therefore, we deduce that 



-w'L- -w' + \w A 



- f 



Cll II i 

W , + 



r 

[A+(f + - f + ) + B(ti - f_)]f + + [A + {w\ - t\) + B(w 2 _ - t 2 _)]w + + 0(r- 

0(r- 8 ), 





-;('♦ 




r 8 / 








"2/+ 



r 6 / 



w , if 1 H — 



0(r" 



-AW+ - BW^_ + A% + B- 



i.e. 



-W? - -W' + \w + + AW + + MW_ = 0(r~ 6 ) for r > R. 



Similarly, we have 
-W" - 



-W'_ + — WL + CW+ + OW_ = 0(r" 6 ) for r > R, 



with A, B, C, D defined as same as in the above process of super-sub solution in different 

cases. Therefore, -(rWL)' = 0(r~ 6 ). On the other hand, notice that for each k > 1 there 

r 

exists a point r k e (k, 2k) such that W' ± (r k ) = (W±(2fc) - W±{k))/k = 0(Ar 7 ) = 0(r^ 7 ), 
then we have that 

\-rW' ± {r) + r k W' ± {r k )\ = 



{rW' ± )'rdr 



<C 2 



r dr < — — 



with constant C2 > 0. Together with r k W' ± (r k ) = 0(r k 6 ) — » as k — >■ 00, we have that 

|rW^_( r )l — ~ f° r r > R sufficiently large, i.e. |M^4( r )l — ~ f° r r — ^ sufficiently large. 

r - C r 

Hence, we deduce that |/4( r ) + — "5" f° r a ^ r > R. □ 



4 Monotonicity 



Next we will present the proof on the monotonicity of the radial solutions. First, we define 
the spaces as in [A BM1] : 

Xo^H^oo^rdr), 



X„ 



00 u 2 

u & X : I —rdr < 00 
r 



u 



x n 



[u'f + u z + — vr 



n 



rdr. 



Of course the spaces X n , n^O, are all equivalent, but we define them this way for notational 
convenience. It is not difficult to show (see |AB] ) that for \n\ > 1, X n is continuously 
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embedded in the space of continuous functions on (0, oo) which vanish at r = and r — > oo, 
and that C£°((0, oo)) is dense in X\. It is possible to define a global variational framework for 
the equivariant problems in affine spaces based on X n+ , X n _ to prove existence of solutions. 



The energy is the same as in (16), except it must be "renormalized" to prevent divergence 



of the -^t~ term at infinity. Here we are only interested in the (formal) second variation of 
this renormalized energy, 



(f+,f-)[u+,u-) :-- 



r I \2 , / / \2 i n + 2 

fu_J +{u_) +~^u + 



n 



rdr 



POO 

/ {[A + (fl - t 2 + ) + B{f_ - tl)]u\ + [A.(f 2 ~ tl) + B{fl ~ t 2 + )]ul}rdr 
2(A + flul + A_f_u 2 _ + 2Bf + f_u + u_)rdr, 



D Er, 



+ 

" /n (63) 
defined for [«+,«_] G X n+ x X n _. 

We have the following fact about radial solutions: 

Lemma 4.1. For any n± 6 Z, if [/+, /_] is the (unique) radial solution of 

D 2 E n+ Mf+J-)[u+,u^] > for all [«+,«_] G X n+ x X n _ \ {[0,0]}. 

In other words, the radial solutions are non-degenerate local minimizers of the renor- 
malized energy. An analogous statement for the Ginzburg-Landau equation with magnetic 
field was derived in |AB GJ. and this observation then became the main step in the proof of 
uniqueness of equivariant solutions proved there. 



Proof. We follow |ABG| . and note that 

/-i 2 



u(r) 

W) 



(u 



A2 



UU 1 f o 



iff 
f 2 



W? - 



r ■ 



Let u± G Co°((0, oo)) (if n± = 0, take u± G C£°([0,oo)) instead). Then [y^, j-] gives an 
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admissible test function in the weak form of ([3]), 



Q = DE n+jn _(f+J-) 







.A 



r u 2 , u 2 1 



/+'/- 



u 



2 ,1 , _<■ U + , a I f2 4.2 - - tt - 



rdr 



+5 



n 2 



./ \2 



( f j 



+^+(/+ - t 2 + )u 2 + + A_(fl - t 2 _)u 2 _ 



n 



n 



H ttU, H -u_ 



+B [Ul ~ t 



(f_-t 2 _)u 2 + ] \rdr. 



2 +)u 2_ , ,,2 ^,.2 



After the regrouping, we get the following useful identity: 



n 2 , n n 2 



( I \2 , / / \2 , '"+ 2 i - 2 
(«+) + («_) +^-« + + ^-«_ 

+ A + (/ 2 - t> 2 + + A_(/ 2 - * 2 > 2 - + 5 [(/ 2 - t> 2 + (f 2 - t 2 _)u 2 + ] } rdr 



-i 2 



/ 2 



-, 2' 



rdr > 0. 



Comparing the formula of the second variation D 2 E n+ ^ n _ (/+, /_) and the above identity, and 



using the fact that A s > is the smallest eigenvalue of the matrix 
that 



A + B 
B A_ 



, we obtain 



fl 



U 



fi 



+ f- 



u 

T- 



2{A + f 2 + u 2 + + A-fiuL + 2Bf + f^u + u_) } rdr 



f2„.2 



/•oo 

> 2 / (A + f 2 + u 2 + + A_/V + 2Bf+f-u+u-)rdr 

Jo 

poo 

> 2\ s / {f + u 2 + + flu 2 _)rdr 

Jo 



>0, 



(64) 



valid for all u± G C£°((0, oo)) (or, u± G Cq°([0, oo)) if the respective n± = 0). The case 
of general u± G X± (or X , in case one of n± = 0) follows by density. It is clear that 
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D 2 E n+ ^ n _(f + , /_) > 0. If it were zero for some [«+,«_], then we would have f + u + = 
—f-U_ = almost everywhere. Since f±(r) > for r > 0, we conclude that D 2 E n+n _ (/+, /_) > 
as claimed. □ 



We are now ready for the proof of Theorem 1.2 from the introduction. 



Proof of Theorem 1.2 . Parts (i) and (ii) follow from the nondegeneracy of D 2 E n+ ^ n _ above; 



part (hi) will be proven in a different manner, using the asymptotics and a compactness 
argument. 

Let u±(r) := f±(r). Differentiating ([3|, we get 

- fl - -Jl + + [Mfl - 4) + B{fi - a]/; 



+ 2A ± f 2 ± f ± + 2Bf ± f T f T - 



i.e. 



< - \u' ± + T ^±u ± + [A ± {f ± - t\) + Btfl - t\)]u± 



2n 2 

+ 2A±f 2 u ± + 2Bf ± f T u T - -±/ ± = 0. (65) 



Now define v± = min{0, u±} < 0,w± = max{0,M-|-} > 0, then u± = v± + w±. For part 



(i), assume B < 0. We multiply the respective equation in (65) by v±, use the facts v + w + = 
and v_w_ = and integrate by parts. By ^ and (]6]), if n± > 1, u±(r) — f'±(r) > in some 
neighborhood r G (0,5). Thus, in case n± > 1, v± is supported away from r = 0. By the 



proof of Theorem 3.2 we may conclude that v± e X n± . Moreover 



/ 1 \ f°° 
v± I u± H — ) rrfr — — I (v'±) 2 rdr, 

o V r / -/II 



with no boundary terms. In case n± = 0, we have u± G X Q by the regularity of solutions, and 
M ±(0) = /±(0) = 0. The integration by parts formula above again holds with no boundary 
terms in this case as well. Therefore, with all above facts, we have the following equations 



o = [ {{v> ± ? + ^4 + \Mfl - 4) + m - t 2 T M 



2n 2 

+2Afc/£4 + 2Bf ± f T u T v ± - ~^f±v± } rdr . 
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Then, use the facts u± = v± + w± and add above two equations together 



o = / { K) 2 + {v'_f + %v\ + n \vi 



+ [A + (fi - ti) + B(f_ - e)y + + [am 2 - ti) + buI - *+)]«- 

+ 2(A + flv 2 + + A_/V ) + 4Bf + f_v + v_ + 2Bf+f4w„v + + w + v_) 

2n 2 . 2n 2 



.y* 3 y 3 



Compare with the formula of the second variation D 2 E n+>n _(f + , /_) in (64), we obtain 

POO 

= D 2 E n+>n _(f+J-)[v+,V-] + 2B / f+f_(w-V+ + w+v-)rdr 

Jo 

POO poo 

+ / t( 1, 4. + v t)Tdr — 2 / — {n 2 ,f + v + + n 2 _f_vJ)rdr. 
Jo r 2 Jo r 

Since n 2 + f + v + , n 2 _f_V-, W-V + and w + v^ are all negative, together with B < 0, we get that 
O^D 2 ^,^ (/+,/_)[«+,«_] 

/*00 poo 

= — / -5-(^4- + v^_)rdr + 2 / — (n 2 _f + v + + n 2 _f^v)rdr 
Jo r l J Q r d 

POO 

— 2B / f + f_(w_v + + w + vJ)rdr 



o 



< 0, 



which is a contradiction unless u±. = 0, i.e. unless f'±{r) > for all r > 0. This proves (i). 

To prove (ii), assume B > and n + > 1 and n_ = 0. This time we multiply the equation 
of u + by v + and the equation of u_ by t>_, and integrate by parts again. Just as in the 
previous case, W- G Xq, and the boundary term in the integration will all vanish. We obtain 
that: 



o = / i «) 2 + + [^(/J - *i) + - e)]t 



ni + 1 



n 2 , 



+2A + f 2 + v 2 + + 2Bf + f_u_v + - 2^/8 r }> /•(/,• 

and 

o = r |k) 2 + + - * a _ ) + B(f + - ti)^ 



n 2 



+2A_fiwt + 2Bf_f + u + w^ - 2— } rdr . 
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As in the first case, we add the above two equations and compare to the formula of the 
second variation D 2 E n+ ^ n {f +) f_) in (64), we get 



< D 2 E n+ , n _(f + J^)[v + ,wJ\ 

COO J 



(vi + wtjrdr - IB 



o 



+ 2 



oo / n 2 



rr 



-f+v + + —f^w_ rdr 



< - 

< 



:{v 2 + + w 2 _)rdr 



which is a contradiction unless t> H 
This verifies statement (if). 



0, w. 



0, i.e. f + {r) > and f_{r) < for all r > 0. 



It remains to prove part (iii) of Theorem |1.2| Denote /±(r; B) the solution of (15) with 
coefficient B. We recall that from ^ of Theorem f±(r; B) ~ r n± for r near 0, that is, 
each component vanishes exactly to degree n± at r = 0, 



/i m) (0;S) = 0, m = 0,...,n ± -l, /^(Ojfl) >0. 



(66) 



When 5 = the system decouples, and each of / 0j ± = /±(-; 0) solves a rescaled equation for 
the standard Ginzburg-Landau vortices, 
1 ^2 



-f" 



- ~fo,+ + -f /o,+ + ^+(/o >+ - n)/o,+ = 0, for r e (0, oo), 



n 

/o,±(0) = 0, / ,±(r) ^ t± as r oo, / ,±(r) > for all r G [0, oo).) 



/o,_ + — /o,- + ^-(/o,- - ^-)/o,- = 0, for r G (0, oo), 



> 



(67) 



For these solutions, it is well-known that (see [HHJ) 

/o )± (r) > for Vr > 0. (68) 

Now choose B > for which the leading-order term in the asymptotic expansion (given 
in (26)) a± = a±{B ) < 0. From the formula (26), it is clear that a± = a±(B) < for all 
B G [0,Z?o]- By the asymptotic estimate (29), there exists R > 0, which may be chosen 
uniformly for B G [0, Bq], such that 



/±(r; B) > for Vr > R,VB e [0, £„] 
here /4( r i ^) depends on both R and parameter B. 



(69) 
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We next claim that /±(r; B) -> /±(r; 0) = / ,±(r) in C fc ([0, #]) as B ->■ 0, for any k > 0. 



Indeed, from Proposition 1.1 iI>b,±{x) = f±(r; B)e in±e is an entire solution to (JTJ in IR 2 . As 
the right-hand sides of Q are L°° bounded, uniformly for B G [0,B ], by standard elliptic 
estimates it follows that the solutions are bounded in C]^ for any a G (0,1). Moreover, 
the right-hand sides of Q are polynomials, and so by a bootstrap argument we obtain the 
stronger conclusion that ipB,± are uniformly bounded in Cf oc for any k > 0. In particular, 
for any sequence B n — > there is a subsequence such that ipB n ,± — > i>o,± in C fc -norm, for 
any k > 2. By the form of the solutions, we have if)o,± = /o,±('")e m±9 , and by passing to the 
limit in the equations we recognize / 0) ± = /±(s 0) as the (unique) solutions to the decoupled 



Ginzburg-Landau system (67) above. By the uniqueness of the limit, we conclude that 
/±(r; B) — > /±(r; 0) = /o,±(r) in C fc ([0, i?]) as 5 — > (that is, along any sequence B n — > 0), 
as claimed. 

We now complete the proof of the statement (iii) by contradiction: suppose there exist 



sequences Bk — > and G (0, oo) such that f±(rk',Bf.) < 0. By (69), we have that 
< rfc < R, so there must exist a subsequence r^. and r G [0, oo) with r^. — >■ r . We claim 
that r = 0. Indeed, if not, by the C 1 ([0,i?]) convergence proven above, 

/«U(n>) = /±(r ;0) = lim f ± (r kj ;B k .) < 0, 



which contradicts (68). Thus r = 0. 

Since f' ± (rkp B k ) = 0, by applying Rolle's theorem (n± — 1) times we may conclude that 
there exist points f kj G (0,r k .) so that f^ n± \f k .) = 0. By the C n± ([0,R}) convergence, we 
then conclude that 

/ ( ; ± ±) (0)= lim/i n±) (f fcj ;5 fcj ) = 0, 



which contradicts (66) for B = 0. Thus we conclude that f'(r; B) > for all r G (0, oo) and 



for all B G [0, B ], as desired. □ 
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